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Abstract 

Prom single vehicle data a number of new empirical results concerning the 
density-dependence of the velocity distribution and its moments as well as the 
characteristics of their temporal fluctuations have been determined. These are 
utilized for the specification of some fundamental relations of traffic flow and 
compared with existing traffic theories. 
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For the prosperity in industrialized countries, efficient traffic systems are indispensable. 
However, due to an overall increase of mobility and transportation during the last years, the 
capacity of the road infrastructure has been reached. Some cities like Los Angeles and San 
Francisco already suffer from daily traffic collapses and their environmental consequences. 
About 20 percent more fuel consumption and air pollution is caused by impeded traffic and 
stop-and-go traffic. 

For the above mentioned reasons, several models for freeway traffic have been proposed, 
microscopic and macroscopic ones (for an overview cf. Ref. [|l||). These are used for devel- 
oping traffic optimization measures like on-ramp control, variable speed limits or re-routing 
systems For such purposes, the best models must be selected and calibrated to empirical 
traffic relations. However, some relations are difficult to obtain, and the lack of available 
empirical data has caused some stagnation in traffic modeling. 

Further advances will require a close interplay between theoretical and empirical inves- 
tigations 0. On the one hand, empirical findings are necessary to test and calibrate the 
various traffic models. On the other hand, some hardly measurable quantities and relations 
can be reconstructed by means of theoretical relations. 

Therefore, a number of fundamental traffic relations will be presented in the following. 
Until now, little is known about the velocity distribution of vehicles, its variance or skewness. 
A similar thing holds for the functional form of the velocity- density relation or the variance- 
density relation at high densities. Empirical results have also been missing for the fluctuation 
characteristics of the density or average velocity. These gaps will be closed in the following. 
Although the data are varying in detail from one freeway stretch to another, the essential 
conclusions are expected to be universal. 

In a recent paper |j| it has been shown that the traffic dynamics on neighboring lanes 
is strongly correlated. Therefore, it is possible to treat the total freeway cross section in 
an overall way. Consequently, we will only discuss the properties of the lane averages of 
macroscopic traffic quantities. The empirical relations have been evaluated from single 
vehicle data of the Dutch two-lane freeway A9 between Haarlem and Amsterdam (for a 
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sketch cf. Fig. 1 in Ref. ||). These data were detected by induction loops at discrete places 
x of the roadway and include the passage times t a (x), velocities v a (x), and lengths l a (%) 
of the single vehicles a. Consequently, it was possible to calculate the number N(x, t) of 
vehicles which passed the cross section at place x during a time interval [t — T/2,t + T/2], 
the traffic flow 

Q(x,t) :=N(x,t)/T, (1) 
and the macroscopic velocity moments 

{vk):= NikT) £ Mx)}\ (2) 

iv ^, oj t _ T / 2 <t a (x)<t+T/2 

Small values of T are connected with large statistical variations of the data, but large values 
can cause biased results for k > 2 |J. Values between 0.5 and 2 minutes seem to be the 
best compromise Jl| . The vehicle densities p(x, t) were calculated via the theoretical flow 
formula 

Q(x,t)=p(x,t)V(x,t). (3) 

Other evaluation methods [[|] are discussed in Ref. 

We start with the discussion of the grouped empirical velocity distribution P(v;x,t) 
which was obtained in the usual way: 

p , + n n{x,vi,t) 

P( ^' t):= ^7T (4) 

Here, n(x,vi,t) denotes the number of vehicles which pass the cross section at x between 
times t — T/2 and t + T/2 with a velocity v G [vi — A/2, Vi + A/2). The class interval length 
was chosen A = 5km/h. 

In theoretical investigations, the velocity distribution P(v; x, t) has mostly been assumed 
to have the Gaussian form 



P G (v;x,t):= . exp n , K \ n . 5 
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Here, V(x,f) := (v) denotes the average velocity and Q(x,t) := ([v — V(x, t)] 2 ) the velocity 
variance. Assumption (|5|) has been made for two reasons: First, it allows to derive approx- 
imate fluid-dynamic traffic equations from a gas-kinetic level of description @-|7]]. Second, 
analytical results for the velocity distribution are not yet available, even for the stationary 
and spatially homogeneous case. Therefore the question is, whether the Gaussian approxi- 
mation is justified or not. Figure [l] gives a positive answer, at least for the average velocity 
distribution at small and medium densities. In particular, bimodal distributions are not 
observed ||. 

An investigation of the temporal evolution of the velocity distribution is difficult due 
to the large statistical fluctuations (which come from the fact that only a few vehicles 
per velocity class pass the observed freeway cross section during the short time period T). 
Therefore, we will study a macroscopic (aggregated) quantity instead, namely the temporal 
variation of the skewness 

_ ([v-V(x,t)f) _ (v*}-3(v)(v 2 } + 2(v) 3 
71 ' J ' [9(x,t)] 3 / 2 [0(M)] 3/2 ' {) 

This can be interpreted as a dimensionless measure of asymmetry (cf. Fig. [2]). Figure |3| 
shows that the skewness mainly varies between —0.5 and 0.5. The deviation from is neither 
systematic nor significant, so that the skewness is normally negligible. This indicates that 
even the time-dependent velocity distribution is approximately Gaussian-shaped ||. 

Now it will be investigated how the average velocity V and the variance G depend on the 
vehicle density p (cf. Figs. ^ and |5]). The problem is that the data for high vehicle densities 
are missing. However, for computer simulations of the traffic dynamics the corresponding 
functional relations need to be specified. This can be done by means of theoretical results. 
For the average velocity and variance on freeways with speed limits, recent gas-kinetic traffic 
models || imply the following implicit equilibrium relations (indicated by a subscript "e"), 
if the skewness is neglected (cf. Fig. [|): 

KM^,- , rMI '-* g ^' (7) 

1 - P/Pmax ~ pT r V e {p) 
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e e (p) = A(p)[v e ( P ) 2 + Q e (p)}, i.e. e e ( P ) = A [ pW f^ ■ (8) 

Herein, Vq denotes the average desired speed (or free speed), r(p) is the effective density- 
dependent relaxation time of acceleration maneuvers, p(p) means the probability of immediate 
overtaking. Moreover, p ma x denotes the maximum vehicle density, T r the reaction time, and 
A(p) with < A(p) <C 1 the relative individual velocity fluctuation during the time interval 
r(p) II- 

According to relation the equilibrium variance vanishes when the average velocity 
becomes zero. This consistency condition is not met by all traffic models (cf. Ref. |fL0|1). 
In addition, we expect that the average velocity vanishes at the maximum vehicle density 
Pmax- Therefore, in the limit p — > p max we must have the proportionality relation 

t(p)[i - P(p)}p ^jfv K 1 - — - p t ^p) > ( ) 

1 ~ ^{P) Pr, 



■'max 



the proportionality factor being Vq. Whereas the overtaking probability p(p) is expected 
to vanish for p —>■ p ma x, the relaxation time r(p) and the fluctuation parameter A(p) are 
assumed to remain finite |pl |. Therefore, the ansatz V e (p) oc (1 — p/p m ax) leads to (3 = 1 



and 



Ve{ P )= P T m ( ax f 2 for p*p max . (10) 

1 r \p r: 



J max , 



This is a very interesting discovery, since many researchers believed that the average velocity 
approaches the p-axis horizontally. In addition, we find that O e (p) oc (1 — p/p max ) 2 for 

P * Pmax- 

Our remaining task is to specify the parameters p m ax and T r . From other measurements 



it is known that p ma x lies between 160 and 180 vehicles per kilometer and lane [0. The 
reaction time T r for expected events is at least 0.7 seconds fl3| . A good fit of the data results 
for 



Pmax = 160 vehicles/km lane, T r = 0.8s 



(11) 
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(cf. Fig. U). In addition we can conclude from (0) that the velocity-density relation V e (p) 
of a multi-lane freeway should start horizontally, since the probability of overtaking p(p) 
should approach the value 1 at very small densities p ~ 0. 

However, it is not only possible to reconstruct the functional forms of the velocity-density 
relation V e (p) and the variance-density relation O e (p)- From these we can also determine 
the dependence of the model functions A(p) and r(p)[l — p(p)} by means of the theoretical 
relations (0) and (J|). The result for the diffusion strength A(p) is depicted in Figure |[ 

Finally, we will investigate the temporal fluctuations of the empirical vehicle density 
p(x, t). Until now, most related studies have been presented theoretical or simulation results. 
It has been claimed that the power spectrum p(x, v) of the density p(x, t) obeys a power law 

z/) oc z/ _l5 , i.e. log p(x, v) = C — 5 log v . (12) 

For 5, the values 1.4 [13], 1.0 [[HJ, or 1.8 [[LlJ] have been found. The empirical results in 



Figure indicate that the exponent S is 2.0 at small frequencies u, otherwise 0.0. Taking 
into account the logarithmic frequency scale, we can conclude that the power spectrum is 
flat for the most part of the frequency range. This corresponds to a white noise. Analogous 
results are found for the power spectrum of the average velocity V(x, t) [l|]. 

In summary, we found that the velocity distribution is approximately Gaussian dis- 
tributed and that its skewness is negligible. We were able to reconstruct the velocity-density 
relation V e (p) and the variance-density relation e (p) by means of theoretical results. This 
allowed the determination of some density-dependent model parameters. The fluctuations 
of the vehicle density could be approximated by a white noise, although a power law with 
exponent 2.0 was found at small frequencies. All these results are necessary for realistic 
traffic simulations. 
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FIG. 1. Comparison of empirical velocity distributions at different densities ( — ) with frequency 
polygons of grouped Gaussian velocity distributions with the same mean value and variance ( — ) . 
A significant deviation of the empirical relations from the respective discrete Gaussian approxi- 
mations is only found at a density of p = 40 vehicles/km lane, where the temporal averages over 
T = 2min may have been too long due to rapid stop-and-go waves H (cf. the mysterious "knee" 
at p ~ 40veh/km in Fig. ||). 




D. Helbing: Fundamentals of traffic Bow, PRE 



10 




-5 -4 -3 -2 -1 1 2 3 4 5 



(v-V) 

Ve 

FIG. 2. Velocity distributions P y (v) := {1 - 7 [3(v - V)/0 1 / 2 - (v - Vf /e 3 / 2 }/6}P G (v) with 
the same average velocity V and variance 0, but different values of the skewness 7 ( — : 7 = 0; 

-: 7 = 1/2; : 7 = 1; • • • : 7 = 2). Obviously, a skewness of | 7 | < 0.5 only leads to minor 

changes compared to the Gaussian distribution ( — ). 
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FIG. 3. Density-dependence of the skewness 7 (• : 1-minute data; O: respective mean values). 
The large variation of the 1-minute data at low densities is due to the small number of vehicles 
which pass a cross section during the time interval T = 1 min, whereas the large variation of their 
mean values at high densities comes from the few 1-minute data, over which could be averaged. 
The 1-minute data of the skewness scatter around the zero line ( — ) and mostly lie between —0.5 
and 0.5. 
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FIG. 4. Relation between average velocity and density (•: 1-minute data; O: respective mean 
values; — : fit function for the equilibrium relation V e (p)). The speed limit is 120km/h ( — ). 
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FIG. 5. Density-dependence of the standard deviation y/Q of the vehicle velocities (•: 1-minute 
data; O: respective mean values; — : fit function for the equilibrium relation ^/6 e (p)). 
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FIG. 6. Density-dependence of the fluctuation strength A(p), which is a measure for the relative 
velocity variation during a time interval r(p). Its maximum at medium densities indicates that 
velocity fluctuations are particularly large in the region of unstable traffic flow. 
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FIG. 7. The power spectrum of the time-dependent vehicle density p(x, t) follows a power law 
with exponent 5 = 2.0 at very small frequencies v, but it is flat over large parts of the frequency 
range, corresponding to a white noise. 



